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Abstract. Gravitino dark matter, together with thermal leptogenesis, implies an upper bound on the 
masses of superparticles. In the case of broken R-parity the constraints from primordial nucleosyn- 
thesis are naturally satisfied and decaying gravitinos lead to characteristic signatures in high energy 
cosmic rays. Electron and positron fluxes from gravitino decays cannot explain both, the PAMELA 
positron fraction and the electron + positron flux recently measured by Fermi LAT. The observed 
fluxes require astrophysical sources. The measured antiproton flux allows for a sizable contribution 
of decaying gravitinos to the gamma-ray spectrum, in particular a line at an energy below 300 GeV. 
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WHY GRAVITINO DARK MATTER? 

An unequivocal prediction of locally supersymmetric extensions of the Standard Model 
is the gravitino, the gauge fermion of supergravity [1]. Its discovery would be of funda- 
mental importance, as the discovery of the W- and Z-bosons for the electroweak theory. 
Depending on the mechanism of supersymmetry breaking, it can be the lightest super- 
particle (LSP), which makes it a natural dark matter candidate. So far, almost nothing is 
known about the gravitino mass. Examples of considered mass ranges are 

• m 3/2 < 1 keV, corresponding to hot dark matter [2], 

• 1 keV < m 3 / 2 < 15 keV, representing warm dark matter [3], 

• 100 keV < m 3/ / 2 < 10 MeV, a range of cold dark matter, favoured by gauge 
mediation and thermal leptogenesis [4, 5], 

• 100 GeV < m 3/ / 2 < 1 TeV, a range of cold dark matter typical for gaugino and 
gravity mediation, suggested in connection with thermal leptogenesis [6]. 

It is very interesting that the various gravitino mass ranges are strongly correlated with 
the masses of other superparticles and with the history of the cosmological evolution. 
This leads to upper bounds on the gluino mass, lower and upper bounds on the next-to- 
lightest superparticle mass (NLSP) and also upper bounds on the reheating temperature 
in the early universe. Moreover, there is a close connection with baryogenesis. 

It would be very exciting to discover a massive gravitino and to establish in this 
way spontaneously broken local supersymmetry as a fundamental, hidden symmetry of 
nature. This appears impossible for the most popular scenario of heavy gravitinos with 
neutralino LSP. On the contrary, it may be possible to discover a gravitino LSP which 
manifests itself as missing energy in NLSP decays [7, 8]. 



GRAVITINO PROBLEM OR VIRTUE? 



In a supersymmetric plasma at high temperature gravitinos are thermally produced, 
mostly by QCD processes. Their number density n 3 / 2 increases linearly with the re- 
heating temperature, 
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where M p and (X3 are the Planck mass and the QCD fine structure constant, respectively. 
The late decay of heavy gravitinos alters the successful BBN prediction, which implies 
upper bounds for the reheating temperature Tr. The most stringent one from hadronic 
decays reads [9] 

T R < ^(1)x 10 5 GeV. (2) 

This is clearly incompatible with thermal leptogenesis [10] which leads to the lower 
bound T R > 10 9 GeV (cf. [11, 12]). 

The conflict between the upper bound from BBN and the lower bound from leptoge- 
nesis on the reheating temperature can be avoided if the gravitino is the LSP [6] . In this 
case the BBN bounds apply to the NLSP which is quasi-stable. One also has to worry 
about gravitino dark matter which may overdose the universe. 

In general, two processes contribute to the generation of gravitino dark matter. In 
the SuperWIMP mechanism [13, 14] gravitinos are produced in WIMP decays. The 
gravitino mass density is then determined by the NLSP density, 

m 3/2 

^3/2 = ^NLSP , (3) 
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which is independent of the reheating temperature Tr. The BBN constraints require, 
however, rather large NLSP masses [15], which makes it difficult to test this mechanism 
at the LHC. The thermal production of gravitinos is dominated by 2 — > 2 QCD scattering 
processes. Solving the Boltzmann equations one obtains [16, 17] 



~ ,2 r,r( T R \ / 100GeV\ fms(n)Y 



where the two-loop running of the gluino mass from the rehating temperature to the elec- 
troweak scale has been taken into account [18]. It is remarkable that the observed dark 
matter density is obtained for typical SUSY breaking parameters, m 3 n ~ 100 GeV, nig ~ 
1 TeV and a reheating temperature characteristic for leptogenesis, T R ~ / 2 Mp ~ 
10 10 GeV [6]. 

The type and masses of NLSP's consistent with leptogenesis and gravitino dark matter 
are strongly restricted by constraints from BBN (cf. [19]), in particular by the catalyzed 
production of 6 Li in the case of a stau NLSP [20]. If the reheating process, which 
leads to the temperature Tr, is taken into account the constraints are sometimes relaxed. 
For instance, in the case of inflaton decays into right-handed neutrinos, § — > N\N\ — > 
(IH)(IH) the reheating temperature required for leptogenesis is smaller by about 

one order of magnitude [21, 22]. If the Boltzmann equation for lepton asymmetry 
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FIGURE 1. Relic gravitino density as function of the gravitino mass; the GUT scale parameter M3 = 
400 GeV corresponds to the on-shell gluino mass nig ~ 1 TeV. The leptogenesis parameters are M\ = 
10 9 GeV and K which lies in the 'strong washout regime'. From [23]. 



and gravitino density are solved simultaneously, a connection between gravitino mass 
and neutrino parameters is obtained directly, without any reference to the reheating 
temperature. As an example (see Figure 1), for M\ = 10 9 GeV, K = 10 (mi = 0.01 eV) 
and itig = 1 TeV one finds m 3/ / 2 — 6 GeV [23]. 

DECAYING GRAVITINO DARK MATTER 

Nucleosynthesis, leptogenesis and gravitino dark matter can all be consistent in the case 
of a small R-parity breaking [24], which leads to the processes — ► TV^, jUV T , Tj, — > b c t, 
... and also to y/3/2 — * yv. Small R-parity breaking couplings can be induced by B-L 
breaking, 

v 2 

A ~ h^®< 10~ 7 , ~ , (5) 

where h^ e ' d ' are the lepton and down-quark Yukawa couplings, respectively. The 'short' 
NLSP lifetimes, e.g., 

typically lead to NLSP decay before BBN. One finds that BBN, thermal leptogenesis 
and gravitino dark matter are consistent for 10~ 14 < A, A' < 10~ 7 and m 3/ / 2 ~5 GeV 
[24]. Characteristic signals at the LHC can be strongly ionising macroscopic charged 
tracks, followed by a muon track, or a jet and missing energy, corresponding to T — > jUV T 
and T — > XV jx, respectively. 
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FIGURE 2. Boundary conditions (B) with stau NLSP. (a) Contours of constant reheating temperature, 
Tr = (1 — 4) x 10 9 GeV, with = £2dm (solid lines); the gray region is excluded by constraints 
from low-energy experiments; thermal leptogenesis is possible in the yellow and orange regions; the 
orange region is favored by the muon g — 2 anomaly at the 2a level, (b) Contours of constant reheating 
temperature in the m stau — m^ji plane; in the dark gray region, the gravitino is not the LSP. From [18]. 
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The gravitino decay ty/3/2 — ► Y v 1S suppressed both by the Planck mass and small R- 
parity breaking couplings, so that the lifetime is much longer than the age of the universe 
[25], 

Decaying dark matter with lifetime i^(10 26 ) s has become very popular after the recent 
observation of the rise in the cosmic-ray positron fraction by the PAMELA collaboration 
[26]. Note that R-parity breaking scenarios with gravitino lifetimes of this order can be 
realized in various ways [27, 28, 29, 30, 31]. 

The consistency of leptogenesis and gravitino dark matter implies important con- 
straints on the superparticle mass spectrum, which depends on the boundary conditions 
of the supersymmetry breaking parameters at the grand unification (GUT) scale. Typical 
examples are 

(A) m = m 1 / 2 , ao = 0; (B) m = 0, m 1 f 2 , ao = , (8) 

with a bino-like neutralino (A), and a right-handed stau (B) as NLSP, respectively. The 
corresponding upper bounds on the gravitino and stau masses are shown in Figure 2 for 
reheating temperatures in the range Tr = (1 — 4) x 10 9 GeV [18]. Relaxing the boundary 
conditions at the GUT scale, gravitino masses up to 1.4 TeV are possible [32]. 
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FIGURE 3. Antiproton flux for mj, n = 200 GeV in the MED set of propagation parameters that saturate 
the antiproton overproduction bound. Dotted line: antiproton flux from gravitino decays, dashed line: 
secondary antiproton flux from spallation in the case of minimal nuclear cross section, solid line: total 
antiproton flux. The gravitino lifetime is t 3/ / 2 = 7 x 10 26 s. From [34]. 



CONSTRAINTS FROM PAMELA AND FERMI 

The R-parity violating gravitino decays i// 3/ / 2 — > yv,hv,Zv 1 W ± l T lead to interest- 
ing cosmic -ray signatures [33] which we shall discuss in the following in a model- 
independent way based on an operator analysis [34]. The mass scales multiplying the 
non-renormalizable operators are inverse powers of the Planck mass Mp and the super- 
symmetry breaking gravitino mass m 3/ / 2 . This assumes for the masses msM of Standard 
Model particles, the gravitino mass and the masses m so f t of other superparticles the hi- 
erarchy m| M -C m 3 j2 <^ m soft- One easily verifies that the dimension-5 and dimension-6 
operators for the R -parity violating couplings of the gravitino are given by 

^eff = jVr V A^ + (ZitfYBpy + bgW^v) o^^Vk } +h-c. , (9) 

where typically £1.2 = &{\/m^/ 2 ). Note, that in general K and the product K^i^ are 
independent parameters which depend on flavour. 

The dimension-5 operator describes the processes i// 3/ / 2 — » hv,Zv, W ± l J? which, after 
fragmentation of Higgs, Z- and W-bosons, yields continuous gamma and antimatter 
spectra, ij/^n — > yX,pX,e + X, which therefore are strongly correlated. It is remarkable 
that the decay ty/ 3/ / 2 — > yv is controlled by the dimension-6 operator. Hence, the intensity 
of a line in the gamma-ray spectrum is not tied to the continuous part of the spectrum. 

The interstellar antiproton flux from gravitino decay suffers from uncertaintes in the 
determination of the physical parameters in the propagation of charged cosmic rays 
in the diffusive halo, leading to uncertainties in the magnitude of fluxes as large as 
two orders of magnitude at the energies relevant for present antiproton experiments. 
The requirement that the total antiproton flux from gravitino decay be consistent with 
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measurements gives a lower bound on the gravitino mass which strongly depends on the 
choice of the halo model. In the following we shall adopt the MED propagation model, 
which provides the best fit to the B/C ratio and measurements of flux ratios of radioactive 
cosmic-ray species [35]. 

A conservative upper bound on the antiproton flux from gravitinos is obtained by 
demanding that the total flux is not larger than the theoretical uncertainty band of the 
MED propagation model. This means that a 'minimal' dark matter lifetime for the MED 
model can be defined by a scenario where the secondary antiproton flux from spallation 
is 25% smaller than the central value, due to a putative overestimation of the nuclear 
cross sections, and the total antiproton flux saturates the upper limit of the uncertainty 
band which stems from astrophysical uncertainties discussed above. 

Using this prescription one finds the lower bound on the gravitino lifetime T™" ~ 

7 x 10 26 s for = 200 GeV. The corresponding antiproton flux from gravitino 
decay, the secondary antiproton flux from spallation and the total antiproton flux are 
shown in Figure 3 together with the experimental measurements by BESS, IMAX and 
WiZard/CAPRICE, and the uncertainty band from the nuclear cross sections in the 
MED propagation model. The minimal lifetime can be compared with the gravitino 
lifetime needed to explain the PAMELA positron fraction excess. 

It is now straightforward to calculate the positron flux at Earth from gravitino decay 
in the MED propagation model [35]. Note that the sensitivity of the positron fraction to 
the propagation model is fairly mild at the energies where the excess is observed, since 
these positrons are produced within a few kiloparsecs from the Earth and barely suffer 
the effects of diffusion. 

To compare the predictions to the PAMELA results, we calculate the positron fraction, 
defined as the flux ratio <£ e + / (4> e + + 4> e ) . For the background fluxes of primary and 
secondary electrons, as well as secondary positrons, we extract the fluxes from "Model 
0" presented by the Fermi LAT collaboration in [36], which fits well the low energy data 
points of the total electron plus positron flux and the positron fraction, and is similar to 
the MED model for energies above a few GeV [37]. Then, the positron fraction reads 

4> D + M (r)+^ g (r) 

pp/j\ _ e+ V I e+ V / (10) 

d>DM (r) + ^kg (r) +4 ,DM (r) +k <j,bkg (r) > 

where k = is the normalization of the astrophysical contribution to the primary 
electron flux, which is chosen to provide a qualitatively good fit to the data. 

For the PAMELA positron excess to be due to gravitino dark matter decay, the 
gravitino mass must be at least 200 GeV. The decay l^ 3 / 2 — > V7 ± £ =F then has a branching 
ratio of ~ 50%, and the hard leptons that are directly produced in these decays may 
account for the rise in the positron fraction if a significant fraction of these leptons has 
electron or muon flavour. 

Consider first the extreme case that the decays occur purely into electron flavour. For 
m 3 /2 = 200 GeV, the PAMELA excess can then be explained for the gravitino lifetime 

t|, 2 (200) ~ 3.2 x 10 26 s, as illustrated by Figure 4. Note that this lifetime is a factor 2 
smaller than the minimum lifetime which was obtained from the antiproton constraint. 
However, the MIN model and other sets of parameters that yield intermediate values 
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FIGURE 4. 'Electron dominance'. Contribution from dark matter decay to the positron fraction (a), 
and the total electron + positron flux (b), compared with data from PAMELA and HEAT, and ATIC, 
Fermi LAT and HESS, respectively; ;« 3 / 2 = 200 GeV, T3/2 = 3.2 x 10 26 s; for W ± P decays pure electron 
flavour is assumed. The "Model 0" [36] background is used, and for comparison with Fermi LAT data 
25% energy resolution is taken into account. From [34]. 



for the antiproton flux can easily be compatible with both the positron fraction and 
the antiproton-to-proton ratio observed by PAMELA. The situation is very similar for 
m 3/2 = 400 and 600 GeV. 

Figure 4 also shows the predicted total electron + positron flux together with the 
results from Fermi LAT [39] and ATIC [40]. Obviously, the "Model 0" cannot account 
for the present data, and the contribution from gravitino decays makes the discrepancy 
even worse. In particular, the data show no spectral feature expected for decaying 
dark matter. On the other hand, gravitino decays may very well be consistent with the 
measured total electron + positron flux once the background is appropriately adjusted. 
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FIGURE 5. 'Flavour democracy'. Contribution from dark matter decay to the positron fraction (a), 
and the total electron + positron flux (b), compared with data from PAMELA and HEAT, and ATIC, 
Fermi LAT and HESS, respectively; m-^ji = 200 GeV, with the minimal lifetime T3/2 — 7 x 10 26 s; for 
W ± P decays democratic flavour dependence is assumed. The "Model 0" [36] background is used, and 
for comparison with Fermi LAT data 25% energy resolution is taken into account. From [34]. 



This is evident from Figure 5 where the contribution from gravitino decays is shown 
in the theoretically well motivated case of flavour democratic decays. The figure also 
illustrates that, depending on the gravitino mass, the dark matter contribution to the 
PAMELA excess can still be significant. 

An obvious possibility is that both, the total electron + positron flux and the positron 
fraction, are dominated by astrophysical sources. For instance, for the gravitino mass 
m 3/2 = 100 GeV one obtains from the antiproton flux constraint the minimal lifetime 
T™"(100) ~ 1 x 10 27 s. The corresponding contribution from gravitino decays to the 
total electron + positron flux and positron fraction turns out to be indeed negligible 
[34]. Nevertheless, as discussed in the next section, the dark matter contribution to the 
gamma-ray flux can still be sizable. 



PREDICTIONS FOR THE GAMMA-RAY SPECTRUM 



The gamma-ray flux from gravitino dark matter decay receives contributions from the 
decay of gravitinos in the Milky Way halo and at cosmological distances, which can be 
calculated in the standard manner. The halo component dominates, leading to a slightly 
anisotropic gamma-ray flux. 

The gravitino decay produces a continuous spectrum of gamma-rays which is deter- 
mined by the fragmentation of the Higgs boson and the weak gauge bosons. In addition, 
there exists a gamma-ray line at the endpoint of the spectrum with an intensity which 
is model-dependent. For our numerical analysis we use a typical branching ratio in this 
channel, 

/ 200 GeV \ 

RR(V/ 3 /2^vr)=0.02 , (11) 

\ m 3/2 J 

for gravitino masses above 100 GeV. In Figure 6 the predicted diffuse gamma-ray flux 
is shown for m 3 / 2 = 100, 200 GeV and the respective lower bounds on the gravitino 
lifetime. These spectra correspond to upper bounds on the signal in gamma-rays that 
can be expected from gravitino dark matter decay. 

For our analysis, two sets of results are used since the status of the extragalactic 
background is currently unclear. For the background obtained by Moskalenko et al., 
the extragalactic component is described by the power law [38] 
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The earlier analysis by Sreekumar et al. led to a less steep background [41], 
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In Figure 6, the slightly anisotropic halo signal has been averaged over the whole 
sky, excluding a band of ±10° around the Galactic disk. For the energy resultion 
o{E)/E = 15% has been used. 




CONCLUSION 

In supersymmetric theories with small R-parity breaking thermally produced gravitinos 
can account for the observed dark matter, consistent with leptogenesis and nucleosynthe- 
sis. Gravitino decays then contribute to antimatter cosmic rays as well as gamma-rays. 
Gravitino masses below 600 GeV are consistent with universal boundary conditions at 
the GUT scale. 

Gravitino decays into Standard Model particles can be studied in a model-independent 
way by means of an operator analysis. For sufficiently large gravitino masses the 
dimension-5 operator dominates. This means that the branching ratios into hv, Zv and 
are fixed. As a consequence, the gamma-ray flux is essentially determined once 
the antiproton flux is known. On the contrary, the positron flux is model-dependent. 
The gamma-ray line is controlled by the dimension-6 operator. Its strength is model- 
dependent and decreases with increasing gravitino mass. 

Electron and positron fluxes from gravitino decays cannot account for both, the 
PAMELA positron fraction and the electron + positron flux measured by Fermi LAT. For 
gravitino dark matter, the observed fluxes require astrophysical sources. However, de- 
pending on the gravitino mass, the dark matter contribution to the electron and positron 
fluxes can be non-negligable. 

Present data on antiproton cosmic -rays allow for a sizable contribution of gravitino 
dark matter to the gamma-ray spectrum, in particular a line at an energy below 300 GeV. 
Non-observation of such a line would place a lower bound on the gravitino lifetime, and 
hence on the strength of R-parity breaking, restricting possible signatures at the LHC. 
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